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We theoretically demonstrate unusual features of superconducting proximity effect in gapped 
graphene which presents a pseudospin symmetry-broken ferromagnet with a net pseudomagnetiza- 
tion. We find that the presence of a band gap makes the Andreev conductance of graphene super- 
conductor/pseudoferromagnet (S/PF) junction to behave similar to that of a graphene ferromagnet- 
superconductor junction. The energy gap An can enhance the pseudospin inverted Andreev con- 
ductance of S/PF junction to reach a limiting maximum value for A n ^> M> which depending on 
the bias voltage can be larger than the value for the corresponding junction with no energy gap. 
We further demonstrate a damped-oscillatory behavior for the local density of states of the PF 
region of S/PF junction and a long-range crossed Andreev reflection process in PF/S/PF structure 
with antiparallel alignment of pseudomagnetizations of PFs, which confirm that, in this respect, the 
gapped normal graphene behaves like a ferromagnetic graphene. 
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I. INTRODUCTION 

Transmission of low energy electrons through a normal- 
metallic-superconducting (N/S) junction is realized via 
Andreev reflection (AR) 1 through which an electron with 
an energy e (relative to the chemical potential fi) and spin 
polarization a, upon hitting the N/S interface, is retro re- 
flected as a hole with the same energy but opposite spin 
direction —a. This peculiar scattering process results in a 
finite conductance of a N/S junction at the bias voltages 
below the superconducting energy gap eV < A5 2 . When 
a ferromagnetic (F) metal is interfaced to a superconduc- 
tor, the exchange splitting of up and down spin subbands 
leads to a suppression of AR at F/S interface. As the 
result, the subgap Andreev conductance decreases with 
increasing the exchange energy h from its value for N/S 
junction and vanishes for a half metallic ferromagnet with 
h = /i, where all carriers have the same spin 3 . The ex- 
change field causes a momentum shift of order 2h/vp be- 
tween Andreev correlated electron-hole in F metal, that 
is responsible for the spatial damped oscillations of the 
induced superconducting correlations in F metalsSHll 

Novel interesting phenomena arise when N/S and F/S 
proximity structures are realized in graphene, the re- 
cently discovered two dimensional solid of carbon atoms 
with honeycomb lattice structure^^. Graphene has a 
zero-gap semiconducting band structure in which the 
charge carriers behave like 2D massless Dirac fermions 
with a pseudo-relativistic chiral property. The car- 
rier type, [electron-like (n) or hole-like (p)] and its 
density can be tuned by means of electrical gate or 
doping of underlying substrate. Recent experimental 
progresses in proximity-inducing superconductivity in 
graphene by fabrication of transparent contacts between 
a graphene monolayer and a superconductor [see for in- 
stance, Refs. [TTHI3] . has provided a unique possibility 
to study relativistic-like superconductivity and proxim- 
ity effect. Peculiarity of AR in graphene N/S junctions 
has been explained by Beenakker, who predicted the pos- 



sibility for a specular AR in undoped normal graphene, 
and its associated anomaly in Andreev c urrent - volt age 
characteristics of a graphene N/S contactP^. In the 
case of a graphene F/S junction, the situation is dra- 
matically different from common F/S junctions. It has 
been shown that for the exchange energies higher than 
the chemical potential h > /i, a peculiar spin-resolved 
Andreev-Klein process at graphene F/S interface can re- 
sult in an enhancement of the subgap Andreev conduc- 
tance by h, up to the point at which the conductance at 
low voltages eV <C A $ is larger than its value for the cor- 
responding N/S structure^"^. Also, the corresponding 
Andreev-Klein bound states in graphene S/F/S structure 
are respons ible f or the long-range Josephson coupling of 
F graphen d 19 * 20 *. Moreover, specific non-local proxim- 
ity effect takes place in graphene-based superconducting 
heterostructures mediating purely by a non-local process 
known as crossed Andreev reflection (CAR) which cre- 
ates a spatially entangled electron-hole pair. While in 
ordinary non-relativistic systems the small value of CAR 
conductance is canceled by the conductance of elastic 
electron co-tunneling (CT) process, it can be enh ance d 
in ballistic graphene N/S/N and F/S/F structureeP 112 ^. 

The above explained unusual properties of AR in 
graphene arise from its electronic structure which is fun- 
damentally different from that of a metal or semicon- 
ductor. Further peculiarity of graphene comes from the 
fact that its electrons, in addition to the regular spin, 
appear to come endowed with the two quantum de- 
grees of freedom, the so called pseudospin and valley. 
The pseudospin represents the sublattice degree of free- 
dom of the graphene's honeycomb structure, and the 
valley defines the corresponding degree of freedom in 
the reciprocal lattice^^ 2 ^. In this paper, we demon- 
strate still another peculiarity of graphene-based super- 
conducting hybrid structures which is resulted from the 
pseudospin degree of freedom of electrons in graphene 
with a (non-superconducting) gap in its electronic spec- 
trum. The effect of the pseudospin and the valley de- 



grees of freedom has already been proven to be drasti- 
cally impor tant in several quantum transport phenomena 
in graphen d 9 * 1Q * 27 ti ^, such that these additional quan- 
tum numbers have been proposed to be used for the 
controlling electronic devices in the same way as the 
electron spin which is used in spintronics and quantum 
computing^"^. Control of the pseudospin of the elec- 
tron has been explained to be more feasible in both mono- 
layer graphene and bilayer graphene^^4 In particu- 
lar, it was shown that the monolayer graphene with a 
gap in its electronic spectrum and an appropriate dop- 
ing presents a pseudospin symmetry-broken ferromagnet, 
with a finite pseudospin magnetization oriented vertically 
to the graphene plane 41 . Here, we find that the pseu- 
dospin degree of freedom has a determining effect on AR 
and the associated proximity effect in hybrid structures 
of gapped normal graphene regions as pseudoferromag- 
nets (PFs) and a superconductor in S/PF and PF/S/PF 
geometries. 

In graphene PF, due to the possibility for a small chem- 
ical potential /i, an electron from the conduction-band 
can be reflected as a hole in the valance-band, depending 
on the electron energy e, /i and the energy gap A at. In 
terms of the carriers pseudospin vector, this corresponds 
to an inversion of the z component of the pseudospin 
vector. Thus, such a peculiar AR is associated with a 
transition from the n-type carriers to the p-type carriers, 
which is called Klein tunnelin^pSEI] [ n analogous to the 
corresponding effect in relativistic quantum mechanics^. 
We obtain that such a Andreev-Klein process at graphene 
S/PF junction can enhance the amplitude of AR and the 
resulting Andreev conductance by A at. In particular, 
we show that depending on the bias voltage the Andreev 
conductance of weekly doped PF (fi <C A at) can be larger 
than its value for the corresponding graphene S/N junc- 
tion. We also demonstrate that depending on the en- 
ergy e of the incident electron, \i and A at, AR can be 
of retro or specular types, respectively, without or with 
the inversion of the z component of the pseudospin vec- 
tor. We further demonstrate that the local density of 
states (DOS) inside PF has a damped-oscillatory behav- 
ior which approves that the energy gap A at in the band 
structure of normal graphene produces an effect similar 
to the exchange field in F graphene. For PF/S/PF junc- 
tion system, we find that the transport is mediated purely 
by CT process in parallel alignment of pseudomagnetiza- 
tions (PMs) and CAR process in antiparallel configura- 
tion, that is accompanied by pseudospin switching effect 
and confirms the crucial rule of the pseudospin in the 
gapped normal graphene and its similarity to the rule of 
spin in a F graphene. 

This paper is organized as follows. In Sec.[llJ we estab- 
lish the theoretical framework which will be used to in- 
vestigate AR in graphene S/PF junction. We present our 
main findings for the Andreev conductance and the prox- 
imity DOS of the S/PF junction, respectively, in Sees. 
Ill and |IV| Section [V] is devoted to the investigation 
of CAR in PF/S/PF junction. Finally, we present the 
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FIG. 1: (Color online) (a) Schematic illustration of the 
graphene S/PF junction, (b) The energy band diagram of 
the S and PF regions to explain the two cases of Andreev re- 
flection at S/PF interface. Left panel: the energy- momentum 
relation of the highly doped S graphene. Middle and right 
panels: the band structure of the n-doped PF region for two 
values of the chemical potential. Middle (Right) panel shows 
that an incident electron from the conduction-band of PF re- 
gion with a subgap energy < e < /n — An {s > [i + A at) is 
reflected as a hole in the conduction- (valance-) band without 
(with) the inversion of the z component of the pseudospin 
at S/PF interface. cr e _ and (Th+ denote the pseudospin vec- 
tors of incident electron and reflected hole. v e and vu denote 
the velocity vectors of the electron and the hole, moving in 
different directions. 



conclusion in Sec. El 



II. MODEL AND BASIC EQUATIONS 

We consider a wide graphene S/PF junction normal 
to x-axis with highly doped superconducting (S) region 
for x < and n-doped pseudoferromagnetic (PF) re- 
gion for x > [see Fig. [lja)]. Pseudoferromagnetic 
region is a monolayer graphene with an energy gap Aat 
in its electronic band structure, which behaves as a pseu- 
dospin symmetry-broken ferromagnet with a perpendic- 
ular to the plane of graphene pseudomagnetization (PM) 
whose direction is switched by altering the type of dop- 
ing between n and p. The magnitude of PM depends 
on the ratio of the chemical potential to the energy gap 
/i/Ajv, such that the vertical pseudomagnetization per 
electron PM Z /N takes its maximum value PM Z /N = 1 
for ji ~ A at, decreases by increasing ji/ An and goes to 
zero for highly doped gapped graphene (/i ^> A at) 41 . The 
most common practical approach to induce a band gap 
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is to realize graphene on top of the appropriat e subs trate 
like SiC which breaks the sublattice symmetry * 43 * 44 !. The 
S part can be produce d by depositing S electrode on top 
of the graphene sheet! 11 " 1 3 1 In this region A n = and the 
superconducting correlations are characterized by the su- 
perconducting pair potential (order parameter) As which 
is taken to be real and constant. To study AR at S/PF 
interface within the scattering formalism, we first con- 
struct the quasiparticle wave functions that participate 
in the scattering processes. In order to describe the su- 
perconducting correlations between relativistic electrons 
and holes of different valleys, we adopt Dirac-Bogoliubov- 
de Gennes (DBdG) equation^ which has the form 



where 



H — n A s 



H 



where 



H = v F (a.p) + A N a z -U(r), 



(1) 



(2) 



is the two-dimensional Dirac Hamiltonian with an energy 
gap, e is the excitation energy and U (r) the electrostatic 
potential is taken to be Uq ^> \i in S region and U = in 
PF region. The electron and the hole wave functions, u 
and v, are two-component spinors of the form (^a^b) 
and cr = (cr^, <j y , <j z ) is the vector of the Pauli matri- 
ces operating in the space of the two sublattices of the 
honeycomb lattice structure. 

An incident electron of the conduction-band from right 
to S/PF interface with a subgap energy e < As can be 
either normally reflected as an electron in the conduction- 
band or Andreev reflected as a hole in the conduction- or 
the valance-band. As illustrated in Fig. []Jb), as long as 
< s < (1 — An the reflected hole is an empty state in the 
conduction-band and AR is retro (middle panel), while 
for e > fx + A n it is an empty state in the valance-band 
and AR is specular, if An < As (right panel). The left 
panel of Fig. [ljb) shows the energy-momentum relation 
of the highly doped S graphene^. The importance of 
AR near the Fermi level imposes the condition of An < 
A s on size of the energy gap An- The retro reflection 
dominates if \i ^> As + A at, while the specular reflection 
dominates if /i <C A^ — An- Using the same method as 
in Ref. [4TJ the pseudospin of the incident electron and 
the reflected hole of the conduction- (valance-) band are 
obtained as, 



(cr(k)) e c = \ (~~^:) 2 (-cosa e x + sina e y) 




(3) 



Aa^ 



(a(k))^ v) = \ (-^ZT^) 2 (-cos a h x±sma h y) 
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(4) 



a e (h) = arcsin ( 



hvq 
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indicates the angle of propagation of the electron (hole), 
at a transverse momentum q with energy-momentum re- 
lation 



C = -M+ \/ A iV 2 + (M&e|) 2 

for the conduction-band electron and 



c(v) 



A N 2 



(hv\k h \ 



(6) 



(7) 



for the hole from the conduction- (valance-) band. Here, 
the superscripts ± denote the two directions of propa- 
gation along the x-axis. As can be seen from the above 
equations when an electron from the conduction-band is 
reflected as a hole in the valance-band, the sign of the 
gap-induced z component of the pseudospin vector (a z ) 
is changed, while in the case of the conduction-band hole, 
it retains its sign. This is shown schematically in the mid- 
dle (right) panel of Fig. [TJb) for the case of Andreev re- 
flected hole from the conduction- (valance-) band without 
(with) the inversion of (a z ) upon AR at S/PF interface. 
Thus, for the incident electron and the reflected hole be- 
ing from different types of bands, we have an inversion 
of the z component of the pseudospin vector upon AR 
at S/PF interface. In the following we will show how 
the pseudospin (<j z ) inversion by AR leads to peculiar 
properties of S/PF and PF/S/PF systems. 

Denoting the amplitudes of normal and Andreev re- 
flections, respectively, by r c ( v \ and va, c (v)i the wave func- 
tions inside PF and S regions are written as 



^ = r c ~ + r c{v) ^ + r At<v) , (8) 
^s=a^ s+ + b^ s ~, (9) 

where ± and ip s± are the solutions of DBdG equa- 
tion for the quasiparticles inside PF and S regions, 
respectively, and the two cases of the Andreev re- 
flected holes from the conduction- (valance-) band with- 
out (with) the inversion of (a z ) are denoted by c(v) in 

ippp . Matching the wave functions of PF and S regions 
at the interface x = 0, we obtain the normal and Andreev 
reflections amplitudes as 

= a c{v) - e ia °+*° + e ia »+*»[l - a <v) e ia °+*°] 

r A , hc{v) 

i A.c(v) — / e 

v ' ycos a e cos ah 



i(a e -a h )/2 ^ _|_ e 2ia h \ 



p (11) 



where 



e{h) 



arcsinh ( 



(fi±s)i-A N 2 



ft = arccos (e/As), 



(13) 
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= (itan/ 



and b c ( v j = sec f3 (— icsc/3). Using 
the above-found reflection amplitudes and wave func- 
tions, we investigate the Andreev conductance of S/PF 
interface and the proximity DOS inside the PF region, 
respectively, in the next two sections. 



III. ANDREEV CONDUCTANCE 

In this section, we evaluate the Andreev conductance 
of S/PF junction by using Blonder-Tinkham-Klapwijk 
(BTK) formula^, 

G c (v) =G (1 - \r c ( v )\ 2 + |r4,c( v )| 2 )cosa e da e , (14) 
Jo 



4e 2 ~ 

Go = -Y^ieV), N(e 



, (15) 



where we put e = eV at zero temperature. The quantity 
Go is the ballistic conductance of N transverse modes in 
a sheet of gapped graphene of width W and 



a c = arcsin (* 



'(/i-s) 2 



(16) 



is the critical angle of incidence above which the Andreev 
reflected waves become evanescent and do not contribute 
to any transport of charge. 

Dependence of the resulting Andreev conductance 
G/Gq on the ratio Ajy/fi is presented in Fig. [2^a) 
for An /As = 0.1 and three different bias voltages 
eV/A s = 0,0.5,1. For A^ < /V(l + eV/A N ), the con- 
ductance decreases monotonically with An //a. In this 
interval, the incident electron and the reflected hole are 
from the conduction-band and therefore AR is without 
the inversion of (a z ) [see Fig. |2jb)]. The density of 
states of the conduction-band hole decreases by increas- 
ing Ajv//^- Thus, the amplitude of AR and hence the An- 
dreev conductance decreases with An / /i and goes to zero 
at An = /i/(l + eV/ An), where the density of states of 
the conduction-band hole vanishes. There is a gap in con- 
ductance for + eV/A N ) < A N < fi/(eV/A N - 1), 
which decreases with eV/As and goes towards smaller 
Ajv//i- For An > ^/(cV/An — 1), the pseudospin (a z ) 
inverted Andreev conductance increases monotonically 
with An /n- In this regime, the transport is between 
the conduction and the valance-band and the incident 
electron of the conduction-band is reflected as a hole in 
the valance band. So the pseudospin (a z ) of the reflected 
hole changes sign [see Fig. [2^b)] and the density of states 
of the hole increases with A^r/V, resulting in an enhanc- 
ing Andreev conductance. In this case, the pseudospin 
(a z ) inverted AR is associated with a Klein tunneling 
of the n-type carriers to the p-type carriers. The en- 
hancing conductance reaches a limiting maximum value 
which depends on the bias voltage. Importantly we see 
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FIG. 2: (Color online) (a) Dependence of the Andreev conduc- 
tance of graphene S/PF contact on the gap Ajv/m ( m units 
of the chemical potential) at three bias voltages eV/As = 
0,0.5,1. (b) The behavior of the z component of the pseu- 
dospin vector (a z ) versus Ajv//x for the incident electron 
from the conduction-band and the reflected hole from the 
conduction (before the gap) or valance-band (after the gap) 
at two values of the bias voltage eV/As — 0.5,1, when 
A N / As = 0.1. 



that depending on the value of the bias voltage, the lim- 
iting value of G/Go for A at ^> M can be larger than its 
value for the corresponding S/N structure (An <C /i). 

For eV/ As = 0, we only have AR without (a z ) inver- 
sion, while for finite bias voltages AR without the inver- 
sion of (a z ) changes to the inverting type when A^/V 
goes across the conductance gap. In order to explain 
the behavior of the carriers pseudospin vector in AR 
process, we have plotted the An / H dependence of the 
z component of the pseudospin vector (a z ) for the inci- 
dent electron e and the reflected hole h in Fig. [2^b) when 
eV/ As = 0.5, 1. In the limit of An <C /i, (a z ) of the elec- 
tron and the hole are zero and the system behaves like a 
graphene S/N structure. Increasing An / H leads to the 
out of plane component for the pseudospin vector of the 
electron and the hole such that (a z ) of the hole increases 
more than that of the electron. In this regime, the con- 
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FIG. 3: (Color online) Plot of the Andreev conductance versus 
Ajv//x for different values of Ajv/As, when eV/As = 1. 



FIG. 4: (Color online) The behavior of the proximity density 
of states (DOS) inside the PF region versus x/Xf for different 
values of Ajv//x, when s/As = 0.5 and A at/ As =0.1. 



ductance decreases monotonically with A at / \i and goes 
to zero at A^ = + eV/ A^)- As is shown in Fig. 
[5|b) , the absence of hole states inside the PF gap causes 
a gap in conductance. For An > /x/(eV/Ajv — 1), (&z) 
of the hole changes sign and decreases with Ajy/fi. In 
this case, the Andreev conductance increases monoton- 
ically from zero and reaches a limiting maximum value 
for A at ^> \i (the limit of specular AR), where (a z ) of the 
electron and the hole have equal magnitudes and differ- 
ent signs. According to Eq. (|3|, the magnitude of (a z ) 
for the electron decreases with eV/As such that in the 
limit of An ^> /i, it reaches 0.2 and 0.1 for eV/ As = 0.5 
and 1, respectively. 

Also the behavior of the Andreev conductance versus 
Ajv//i is shown in Fig. [3] for different values of A n /A s, 
when eV/As = 1. It is seen that the limiting value 
of G/Go increases by increasing A at /As for An <C \i 
and tends to the corresponding value of a retro type AR 
(G/Go = 2) as A N -> As, while for A N > [i it de- 
creases from its value for a specular AR in corresponding 
S/N structure An <C As (G/Go = 4/3) and vanishes for 
An > A^. Also the Andreev conductance gap is getting 
broadened by increasing A n /As . 

So the behavior of the Andreev conductance with 
Ajv//i is similar to that of a graphene F/S junction with 
where AR of n-n type carriers for < \i changes 
to the Andreev-Klein reflection of the n-p type carriers 
for h > /x 16 . This tells us that the energy gap A at in 
the band structure of normal graphene behaves like an 
exchange energy in F graphene and enhances the sub- 
gap Andreev conductance of S/PF junction, which is ac- 
companied by the inversion of the z component of the 
pseudospin vector for the reflected hole relative to the 
incident electron. 



IV. LOCAL DENSITY OF STATES 

Let us now study the proximity effect in S/PF junc- 
tion by focusing on an experimentally accessible quantity, 
the local density of states (DOS) inside the PF region. 
To find the proximity DOS as a function of energy and 
position, we use the formula^ 

N(e,r) = J2\Mr)\ 2 5(e(k)-e), (17) 

k 

where ^(r) corresponds to the eigenfunction of energy 
e(k) and the sum is over all states with the wave vectors 
k. Replacing the wave function of Eq. ([8| in the above 
equation, we find the total subgap DOS inside the PF 
region for two cases of AR without (with) (a z ) inversion 
as, 

N« v \e,x) 1 / A N 
N (e) 4y V + ^ j 

/tt/2 
|^(r)p ( ^| 2 cos 2 a e da e , (18) 
-tt/2 



^ UV c; — V -L^y 

(^ F )y(/i + £) 2 -A A r 2 

Here, N (e) is the DOS of a PF layer. Fig. [I] shows the 
behavior of the proximity DOS inside the PF region in 
terms of the dimensionless distance x/Xp (Xf = hvp/fj) 
for different values of Aat/m, when e/As = 0.5 and 
Aat/As = 0.1. The results are for the case of AR with 
(a z ) inversion at S/PF interface. It is seen that N(x) 
decays rather quickly close to the interface with a slope 
which increases by increasing An / '/i- We can see that 
the length scale characterizing this fast decay decreases 



6 



by increasing An //a. Slightly away from the interface, a 
much slower oscillating behavior can be seen with the 
wavelength which clearly decreases with Aj^/fi. Also 
the length scale which characterizes these oscillations de- 
creases with Ajsf/fi. 

So we conclude that there are two phenomena to con- 
sider in describing the spatial variations of N(x), when 
an energy gap is present in the band structure of nor- 
mal graphene. The first is the short distance decay at 
the interface and the other important phenomenon is the 
damped oscillation of N(x), caused by the momentum 
shift between Andreev correlated electron-hole with op- 
posite (<j z ) directions. As can be seen from Fig. [4j the 
period of oscillations is determined by hyp /An, which 
is similar to a S/F structure where the period of DOS 
oscillations in the ballistic limit is given hvp/JW^Q. This 
shows the similarity of the effect of an spin-splitting ex- 
change field h with the energy gap An, which behaves 
as a pseudospin-splitting field [see Eq. The gen- 

eral method to detect these oscillations is tunneling spec- 
troscopy, which is widely used to probe the DOS oscilla- 
tions in S/F structures with different thicknesses of the 
F layer [see e.g. Ref. 4 . In our case, the local scan- 
ning of the surface of the PF region by scanning tun- 
neling microscopy (STM) is the appropriate method to 
probe the DOS oscillations. We note that the imper- 
fectness of the surface at S/PF junction cannot change 
the qualitative behavior of the results. It may lead to 
weakening the proximity effect and changing the results 
quantitatively 5 46 . So the spatially-damped oscillatory 
behavior of the DOS inside the PF region confirms that 
the energy gap An in the band structure of normal 
graphene produces an effect similar to the exchange field 
in F graphene. 
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FIG. 5: (Color online) Plots of the conductance of CAR 
process in antiparallel configuration of pseudomagnetizations 
versus the length of the S region for different two values of 
An/ As, when fx/ An — 1-1 and eV = ji — An- 



V. CROSSED ANDREEV REFLECTION IN 
PF/S/PF STRUCTURE 

We further study the non-local quantum transport and 
CAR in PF/S/PF junction which constitutes a supercon- 
ducting pseudospin valve structure. In CAR process an 
electron excitation and a hole excitation at two separate 
PF leads are coupled by means of Andreev scattering pro- 
cesses at two spatially distinct interfaces. The supercon- 
ducting pseudospin valve consists of two PF regions with 
a tunable direction of PM, which are connected through 
a S region of length L. The configuration of PMs in the 
pseudospin valve can be changed from parallel (P) to an- 
tiparallel (AP) by fixing the type of doping of one region 
and changing the type of the doping in the other region. 
Using the solutions of DBdG equation for the quasipar- 
ticles of PF and S regions, we write the wave functions 
inside the two PF and S regions of P and AP configu- 
rations, within the scattering formalism. Matching the 
wave functions at the two interfaces, we calculate the 
normal and Andreev reflection amplitudes in the left PF 
region and the transmission amplitudes of the electron 



and the hole into the right PF region of both P and AP 
configurations. Replacing the reflection and transmission 
amplitudes in BTK formula, we obtain the conductance 
of AR, CT, and CAR processes for P and AP alignments 
of PMs. We find that for all incoming waves with two 
bias voltages eV = ±(/z — A at), AR process is suppressed 
and the cross-conductance in the right PF region depends 
crucially on the configuration of PMs in the two PF re- 
gions. We find that the transport is mediated purely 
by CT in P configuration and changes to the pure CAR 
in the low energy regime, by reversing the direction of 
PM in the right PF region. This suggests a pseudospin 
switching effect between the pure CT and pure CAR in 
PF/S/PF structure, which can be seen from Eq. ^ ( 
Eq. (Hi) for the right going conduction- (valance-) band 
electron (hole) of n- (p-) doped PF region by replacing 
— cosa e with cosa e (/i with — fi). 

Fig. [5] shows the behavior of the conductance of CAR 
process in AP configuration versus the length of the S 
region for different values of An/ As, when /i/ 'An = 1.1 
and eV = \i — An- It is seen that the CAR conduc- 
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tance has an oscillatory behavior with L/£ and increases 
by increasing An/ As from its value for the correspond- 
ing graphene N/S/N structure. Also we can see that in 
contrast to the graphene N/S/N structure, CAR process 
is present for long lengths of the S region. This effect is 
similar to graphene F/S/F structure^ and approves that 
the gapped normal graphene behaves like a F graphene. 

VI. CONCLUSION 

In conclusion, we have investigated proximity effect 
in graphene-based hybrid structures of superconductors 
and gapped regions as pseudoferromagnets. A gapped 
graphene is in a sublattice pseudospin symmetry-broken 
state with a net pseudomagnetization oriented perpen- 
dicularly to the plane of graphene. We have found that 
upon a certain condition, Andreev reflection of an elec- 
tron from a S/PF interface is associated with an inver- 
sion of the perpendicular component of its pseudospin, 
and that this has important consequences for the prox- 
imity effect. For a S/PF junction, we have found that, 
the Andreev-Klein reflection can enhance the pseudospin 
inverted Andreev conductance by the energy gap A at to 
reach a limiting maximum value for An ^> /i, which de- 
pends on the bias voltage and can be larger than the 
value for the corresponding junction with no energy gap 
(An <C /i). This is similar to the behavior of Andreev 



conductance with the exchange energy h in a graphene 
ferromagnet-superconductor junction. We have further 
studied the proximity density of states (DOS) in PF side 
of S/PF contact, which exhibit a damped-oscillatory be- 
havior as a function of the distance from the interface. 
The period of DOS oscillations is found to be inversely 
proportional to the energy gap An- The proximity DOS 
in ferromagnetic graphene shows similar spatial oscilla- 
tions with a period determined by 1/h. For a PF/S/PF 
structure, we have found a pseudospin switching effect in 
which the subgap transport of electrons can be switched 
from a purely elastic electron co-tunneling process to the 
pure crossed Andreev reflection by changing the align- 
ment of the pseudomagnetizations of PF regions from 
parallel to antiparallel configurations. This is again simi- 
lar to the behavior of the corresponding superconducting 
structure with ferromagnetic graphene. This confirms 
that, in this respect, the effect of the sublattice pseu- 
dospin degree of freedom in gapped graphene is as im- 
portant as the spin in a ferromagnetic graphene. 
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